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ON  TUB  STRUCTURE  OF  SHOCK  WAVES 


by 

K.  Khandrikh 

Abstract 

The  structure  of  shock  waves  in  a  thermally  conducting  viscous 
gas  is  examined.  A  "law  of  conservation  of  irreversible  energy  flow" 
within  a  shock  front  is  derived.  Approximate  analytic  expressions 
for  the  integral  curve  are  obtained  for  the  cases  where  Prandtl's  num¬ 
ber  P  1  and  P<1  and  for  the  velocity  profile  in  a  plasma  with  equal 
electron  and  ion  temperatures  in  the  approximation  P  =  0.  — 

Introduction 


There  are  numerous  methods  for  determining  the  structure  of  a 
shock  wave.  However,  these  numerous  results  are  poorly  suited  to  an 
investigation  of  its  general  properties,  since  the  structure  depends 
ori  several  parameters  (for  example,  the  Mach  number  (M),  the  ratio  of 
specific  heats  (y),  Prandtl's  number  (P)). 

Certain  general  conclusions  about  the  magnitude  of  the  irrever¬ 
sible  flows  and  the  profile  of  a  shock  wave  can  be  obtained  with  the 
help  of  the  "law  of  conservation  of  irreversible  energy  flow,"  which 
is  presented  in  the  first  part  of  this  paper.  One  can  clearly  indi¬ 
cate  in  the  "velocity-temperature"  plane  what  contribution  anisotropic 
pressure  and  heat  flow  make  to  the  irreversible  energy  flow. 

In  the  second  part  of  the  paper  approximate  analytic  expressions 
are  derived  for  the  integral  curve  in  two  cases,  PJ>1  and  P<fl.  •  One 
can  ascertain  by  means  of  those  expressions  the  accuracy  with  which 
the  well-known  integral  curves  for  P=  ra  ,  o  apply  in  the  case  of'  large 
and  small  values  of  P.  They  also  permit  the  refinement  of  the  analytic 
expressions  for  the  integral  curve  in  these  cases.  Small  Prandtl  num¬ 
bers  (  ~0.05)  occur,  for  example,  in  a  plasma  with  equal  electron  and 
ion  temperatures  /2/,  largo  values — (~800),  for  example,  in  neutral 
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carbon  dioxide  gas  /'}/. 

In  the  last  section  an  approximate  analytic  expression  for  the 
velocity  profile  is  discussed  under  the  assumption  of  equality  of  the 
electron  and  ion  temperatures  and  with  the  use  of  the  integral  curve 
for  P  ■=  0. 


Law  of  Conservation  of  irreversible  Energy  Flow 
in  a  Shock  Front 


A  shock  wave,  as  is  well  known,  is  formed  as  a  result  of  nonlinear 
effects,  which  cause  an  increase  in  the  gradients.  The  size  of  these 
nonlinear  effects  depends  only  upon  parameters  characterizing  the  equi¬ 
librium  of  the  gas  and  the  flow  velocity.  In  stationary  shock  waves 
these  nonlinear  effects,  which  cause  an  increase  in  the  gradients,  and 
the  irreversible  processes,  which  decrease  the  gradients,  arc  balanced 
out. 

From  the  equations  of  continuity,  conservation  of  momentum  and 
energy  for  a  stationary  plane  shock  wave  in  an  ideal  gas  (see,  for  ex-' 
ample,  /4/) ,  one  obtains  after  simple  transformations 


V1  ~  1 1  (i) 
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(((denotes  a  quantity  in  front  of  the  shock  wave), /^- is  the  anisotropic 

pressure,  £  the  irreversible  heat  flux,  ,  the  Mach  number, 

v  cp  .  \  W  u  ) 

y  e=  — _  ,  the  ratio  of  specific  heats,  and  the  remaining  symbols  are 

the  usual  ones. 

If  (2)  is  subtracted  from  (3),  then  we  get 

vpxx  —  (Y—  U<7  =  /u(t0~O.  (5) 

With  the  help  of  the  expression  (4)  this  equation  can  be  rewritten  in 
the  form 

u  —  1  ci  *—  1  iif  r:  f ..  ■  i  ivi  ~  1  1  —  — -  —  1  1  1  —  7  1  7 

(6) 
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and 


v-i-i 


(jS  denotes  a  quantity  behind  the  front). 

Equation  (5)  or  (6)  represents  the  desired  "law  of  conservation 
for  irreversible  processes."  The  moduli  of  separate  irreversible  energy 
flows  contribute  to  the  expression  U /JVV — (y — i )<]  i  which  is  responsible 
for  the  weakening  of  the  gradients,  since  it  is  always  the  case  that 
0  and  q-0.  On  the  right-hand  side  of  (5)  or  (6)  there  appears 
the  flux  caused  by  the  nonlinear  effects.  It  is  easy  to  verify  that 
in  the  linear  approximation  in  u  this  flux  is  equal  to  zero. 

In  the  stationary  case  the  moduli  of  both  energy  flows  are  equal 
foi  any  value  of  the  velocity  u .  From  (2),  (3)»  and  (5)  a  simple  clear 
representation  in  the't-u  plane  (Fig.  1)  of  the  irreversible  energy 
flows  and  the  sums  of  their  moduli  follows.  According  to  Fig.  1,  the 
integral 


curve? =  't  (u)  has  no  influence  on  the  "sum"  of  the  irreversible  energy 
flows.  It  only  determines  what  contribution  the  anisotropic  pressure 
fix  and  the  heat  flux  ^  make. 

In  order  to  obtain  the  velocity  and  temperature  profiles,  one  must 
still  know  the  explicit  dependence  of  the  transport  fluxes  on  the  vel¬ 
ocity,  the  temperature,  and  their  gradients.  This  dependence  has  the 
simplest  form  in  the  Navier-Stokes  approximation 

*.---5-’  •‘(n-2-.f— <7> 

If  the  law  of  conservation  of  irreversible  energy  flows  (6)  is  applied, 
then  in  particular  the  following  well-known  results  are  easily  obtained 
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with  the  aid  of  the  Navier-Stokcs  approximation.  In  the  event  of  an 
increase  in  the  transport  coefficients^/^  and  yt  ,  the  gradients  -j—  and 
^  must  decrease |  and  consequently  tho  width  of  the  front  is  increased. 
Since  the  conservation  law  (6)  is  valid  at  each  point  of  the  shock  front, 
this  conclusion  is  true  for  any  segment  of  the  front.  Therefore,  it 
follows  that  an  increase  of  the  transport  coefficients  with  temperature 
causes  the  gradients  to  decrease  more  and  more  in  the  direction  of  the 
end  of  tho  shock  front  and  the  velocity  and  temperature  profiles  to  be 
extended.  In  particular,  if  the  viscosity/*  and  the  thermal  conducti¬ 
vity  “)i.  depend  upon  the  temperature  as  Ts,  then  one  can  obtain  the  de¬ 
sired  profile  from  the  profile  corresponding  to  j  »  O  by  means  of  the 
transformation  dx-wdx'  dx. 

A  similar  law  of  conservation  of  irreversible  energy  flows  can 
also  be  formulated  for  shock  waves  in  which  other  irreversible  pro¬ 
cesses,  such  as  radiative  exchange,  diffusion,  etc.,  play  a  significant 
role.  For  example,  for  a  completely  ionized  plasma  consisting  of  ions 
and  electrons  which  are  moving  with  the  same  velocity,  we  obtain 

v(p<„  +  rj-iy  -')(<?'+  <?')  -  r  ~ IT  MS(‘  ”  y (l  “ z)  *• 

where  <?'■  qr  are  the  corresponding  pressure  and  thermal  flux 

of  the  ions  and  the  electrons.  This  equation  is  also  valid  in  the 
case  where  the  ion  and  electron  temperatures  differ. 

Dividing  (3)  by  (2)  we  get  the  equation  for  the  integral  curve 

?  *=?(«) 

3  1  1  d  ?-  >  ^  u  (8) 

X  "  P  M2’  du  =  >-  ?0 

where  P  cp  is  Prandtl’s  number*  and  cp  is  the  specific  heat  capacity 
at  constant  pressure. 

♦According  to  /5 /  one  may  neglect  the  dependence  of  P  on  tho  tempera¬ 
ture  for  practically  all  possible  laws  of  interaction  between  particles 

J  if  an  ideal  gas  is  considered  (cn  °  constant).  - - 

As  is  well  known, there  are  three  cases  in  which  a  simple  analytic 
expression  for  the  integral  curve?*  =  ?"(u)  is  obtained,  namely  for 
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P  =  eo,  0,  and  3/4  (see,  for  example,  /3/). 

It  is  reasonable  to  suppose  that  a  small  change  in  the  Prandtl 
number  will  result  in  only  an  insignificant  change  in  the  integral 
curve.  Near  the  points  u  ,  u  one  can  easily  prove  the  validity  of 
this  supposition  with  the  aid  of  the  linear  approximation.  On  the  basis 
of  these  considerations  one  can  attempt  to  find  approximate  analytic 
expression  for  the  integral  curve,  if  P  is  insignificantly  different 
from  the  three  cases  mentioned  above,  using  the  known  solutions  as 
the  zeroth  approximations.  In  the  paper  by  Grad  /l/  an  approximate 
method  was  worked  out  for  P  near  3/4.  In  the  present  section  ex¬ 
pressions  are  derived  for  the  integral  curve  in  the  case  of  large 
and  small  P. 


T  —  T»  = 


(9) 


P»l.  '  If  we  write  (8)  in  the  form 
_3  |  t  —  tn  rfr 

4  P  M  ‘  u  ~0u  ' 

then  it  is  obvious  that  for  sufficiently  large  Prandtl  numbers  P  one 

can  neglect  the  right-hand  side  and  thus  obtain  the  zeroth  approxima¬ 


tion'  to  the  integral  curve 


T«»  =  , 


(10) 


Inserting  the  zeroth  approximation  (10)  into  the  right-hand  side  of 
(9) i  we  get  the  first  approximation 


T(l)  LJ I  -  T|0)  —  T|,  rft|0)  '  (11) 

4  P  M1  u  dn  ' 

and  analogously 

T<n>  _  _  3  1  I  t1"-11  —T0  dt*'1-1* 

“  4  P  At*  u  ~  du  '  (12* 

It  is  necessary  for  the  convergence  of  the  method  that  all  thejh^L-for 

Ou  .  . 

ii  in  the  segment  |Woj  remain  finite.  With  the  help  of  equations  ) 
and  (l0)-(12)t  this  condition  is  fulfilled  for  any  values  of  u,»  Mf  and 
Y  •  In  the  first  approximation  to  the  integral  curve  one  obtains  in 
an  obvious  fashion,  for  example,  the  following  expression: 


TO)  - 1-  =  -f  J-  —■  (V  - 1)  -f--  II  +  yAI*  (1  -  «))  X 
X  £  yAJ5  —  — — -  \M3  (1  —  ii)  —  yJ  - 


P«l,  In  this  case  we  write  equation  (8)  in  the  form 
f—  T0  ~~PMht(x— 1«)  -j- 

3  ax 

du 
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and  analogously  obtain 


,  *<0'  f=  T„ 

t<"  - 10  4  PM'uUm  -  xv )  — !_ 

3  *(”> 


du 
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T<«>  -  T„  =:  -i  pM'-u  (T(»-  I)  _  T  )  . 

3  *<»-.!>  ' 

dti 

dx{n) 

In  order  of  this  method  to  converge  it. is  necessary  that  all  the  ' 
differ  fonn  zero.  This  condition  is  not  always  satisfied.  It  is  well 
known  (see,  for  example,  /4/)  that  under  the  condition 

_L  3y_-^i</M, 


V  3-y 
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on  the  interval  ([/'((,  nowhere 

«»  = 


5T(»-'-v) 
r(2-^) 


Y  + 
_ l_ 

Y  + 


there  exists  for  P  =  0  the  so-called  isothermal  jump  and  ((/„,  Up],  There¬ 
fore  we  apply  the  method  discussed  to  the  entire  range  of  velocities 
only  under  the  condition 

JW’C-Livni  (i4) 

Y  3  —  y 

However  one  can  assume  that  in  the  opposite  case  (under  condition  (13) ) 
the  approximations  are  valid  on  the  half-open  segment f[/a ,  i/p].  In  the 
first  approximation  to  the  integral  curve  we  get 
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.  '  3 
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The  detailed  analysis  of  the  convergence  of  the  methods  discussed  is 
fraught  with  well-known  difficulties.  The  region  and  rapidity  of  con¬ 
vergence  depends  upon  all  the  parameters  u,  M,  andy. 

In  Fig.  2  the  first  few  approximations  for  P  =  10;  1 /20 ; y  =  5/3; 
M  ~  6/5  are  shown  as  an  example. 


Velocity  Profile  in  a  Plasma 

In  this  section  we  try  to  obtain  an  approximate  analytic  expres¬ 
sion  for  the  velocity  profile  in  a  plasma  under  certain  simplifying 
assumptions.  According  to  /6/  the  electron  and  ion  temperatures  in  a 
shock  wave  do  not  differ  very  much  from  one  another  for  M<-2.*  There¬ 
fore,  in  what  follows  we  assume  that  the  electron  and  ion  temperauures 
are  equal.  If  the  relative  error  of  the  integral  curve  is  determined 
with  the  help  of  the  expression 

|T»)_t(0)| 

T„-l  ’ 

then  in  the  case  of  a  plasma  with  equal  electron  and  ion  temperatures 
(pps’0.05,  /2/)  we  obtain  for  example,  /l-I„„s«3%  (see  Fig.  2)  when 
M  =  6/5  and  when  M®»1.  This  means  that  'fc'c  -^represents 

a  sufficiently  good  approximation  to  the  integral  curve  for  a  plasma. 


•In  this  case  one  can  also  neglect  the  influence  of  a  balancing  of  the 
electron  and  ion  temperatures  on  ihe  structure  of  the  shock  wave. 


Substituting  the  Navier-Stokes  approximation  (7)  into  (2)  for  £ 

and?-  =  7-  (3)  for  the  integral  curve  and  noting  that-  _fh_  _  ~  tl"~  'and 
%  if.r  tin  <tv 

x  — p~r'2'  *  we  obtain 


a 

«  [ Al-  (I  •—  21)  -h  1  ]  [ ~  Al*  {I  —  I1)  +  <]  *  (15) 
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where 


l 

P 


«„  ys&rm 


(P is  the  mean  free  path  for  momentum  transport. )  It  is  found  that 
one  can  calculate  the  integral  in  (15)  exactly.  The  simple  but  tedious 
calculations  result  in  an  expression  which  we  present  in  an  already 
simplifiei  form,  which  represents  a  good  approximation  for  2, 


j  -  V^'j) 

-  /  "a 


» -I 


4  ,  5  ,  ,  S  ,  0 

T  +  T*,  +  T,,  +  T 


■i - —  arcctg  sj  +  c  («,), 


(16) 


where 


a*  —  I 
3a1  —  5 


Z?ra 
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2'i  3as 


(17) 


The  constant  of  integration  0(uo)  is  chosen  such  that  u  =  uQ  at  £  =  tD, 

where  uD  and  are  arbitrary  quantities.  For  weak  shock  waves  (M-'-l) 

.the  first  two  terms  on  the  right-hand  side  of  (16)  play  the  main  role, 

since  A~B'~  — — ;  .  In  shock  waves  of  arbitrary  strength,  the  first 

alt— I 

term  makes  the  main  contribution  when  u-tuftand  the  second  term  when 
u-*-U|j .  The  third  term  is  important  only  for  the  center  of  the  shock 
wave.  Thus  it  is  evident  that  with  an  increase  in  M  the  profile  is 
extended  to  the  end  of  the  shock  wave  (u  u),  since  the  second  and 
third  terms  in  (16)~M  ,  whereas  the  first  term  is  independent  of  M 
in  the  first  approximation. 

We  recall  once  more  that  the  expression  (16)  gives  the  correct 

profile  in  the  entire  range  «(,]  only  under  condition  (14).  In  the 

opposite  case  (16)  is  applicable  only  in  the  range  [»„,  i/'p)  and  the 

change  ur»u'  occurs  "at  a  single  instant." 

f  f 
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